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Abstract 



A Langevin equation with multiplicative noise is an equation schematically 
of the form dq/dt = — F(q) + e(q)£, where e(q) £ is Gaussian white noise 
whose amplitude e(q) depends on q itself. Such equations are ambiguous, 
and depend on the details of one's convention for discretizing time when 
solving them. I show that these ambiguities are uniquely resolved if the 
system has a known equilibrium distribution exp[— V(q)/T] and if, at some 
more fundamental level, the physics of the system is reversible. I also discuss 
a simple example where this happens, which is the small frequency limit 
of Newton's equation q + e 2 (q)q = -VF(q) + e -1 (q)£ with noise and a 
q-dependent damping term. The resolution does not correspond to simply 
interpreting naive continuum equations in a standard convention, such as 
Stratanovich or Ito. 

I. INTRODUCTION 

A Langevin equation of the form 

<7q=-W(q)+C, (la) 



(Q(t)Q(t'))=2aT6 lJ 6(t-t'), (lb) 

can appear as the effective description of highly overdamped motion of some coordinates q 
in a potential V in contact with a thermal bath. The damping, and the random force term, 
arise from interactions with the thermal bath. By overdamped, I mean that at a slightly 
more fundamental level, the equation of motion might be, for instance, 

q+<rq=-W(q) + C, (2) 

with £ as before. Here, a is a damping coefficient, and £ is a Gaussian thermal noise term. 
In the limit that a is large compared to the inverse time scales of interest, one can ignore 
the q term and so obtain QTa]) . 

In this article, I want to consider the case where (a) the damping coefficient a depends 
on q itself, 



^i(q) Qi = -ViV(q) + (i, (3a) 

(Ci(*)0(O> = 2^(q)T(5(t-O, (3b) 

and (b) it is already know through other means that V(q) is an effective potential that gives 
the equilibrium distribution of q as 

P eq (q) = exp[-V(q)/T] (4) 

in whatever approximation may be relevant to the problem of interest. Equation ([3]) is a 
special case of what are known as Langevin equations with multiplicative noise. In general, 
such equations are notorious for being ambiguous — they are sensitive to exactly how one 
discretizes time. The purpose of this paper is to show that, with a few very general as- 
sumptions, the property (|j) is sufficient to uniquely resolve this ambiguity. Moreover, this 
resolution does not correspond to any standard interpretation (Ito or Stratanovich) of the 
Langevin equation (§). 

Though my arguments will be somewhat general, it helps to have in mind a concrete 
example of an wn-ambiguous description for which @ will be a limiting case. A simple 
example is the analog of (0) with q-dependent damping: 

Qi + o-y(q) 4j = -ViV(q) + Ci, (5) 

with the noise again (|3"15|). As I discuss in Appendix [A|, this equation is unambiguous and 
has equilibrium distribution (Q). 

I can't resist mentioning the application of interest to me personally, which motivated 
this work electroweak baryon number violation in the early universe. Its study requires 
understanding the effective dynamics of fluctuations in weakly-coupled high-temperature 
non-Abelian gauge theories, where there is an effective theory of the form (|3]) at the relevant 
distance and time scales, but where it is much more straightforward to analyze static issues, 
such as (01), than the subtleties of dynamical ones. 



II. A FIRST PASS 

Let's rewrite (§) in the equivalent form 

q=-F(q)+e(q)&, (6a) 

= 2T <%<*(*-*')■ (6b) 

where the q dependence has been scaled out of the noise by defining £ = e(q) £, and the 
matrix e and vector F are 

e(q) = [a(q)]- 1 ^, (7) 
F(q)=a- 1 (q)Vy(q). (8) 
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Two standard conventions for discretizing time, and so removing the ambiguities inherent 
in equations like (H), are the Ito convention^] 

q* - qt-At = -At F It6 (q t _ At ) + e(^ At ) £ t , (9) 

and the Stratanovich convention, 

q* - qt-At = -At F strat (q) + e(q) £ t , (10) 

qt + q*-A* 

q " 2 ' 

where in both cases the discretized noise correlation is 

(i i £ j1f ) = 2TAt6 ij 6 te . (12) 

The specification of F Strat (q) as opposed to F Strat (q(_At) in (|10D is actually irrelevant: it is 
only the q used to evaluate e that causes the difference between these two conventions in 
the At -»• limit. 

If one simply tries setting F It5 = F or F strat = F, the Ito and Stratanovich equations will 
give rise to different physics, and in particular different equilibrium distributions. The two 
conventions are identical in the At — ► limit only if one sets 

Ft 6 = Ft^-Te iad e ja . (13) 

(Here and throughout, I adopt the notation that indices after a comma represent derivatives: 
F itj = dFi/dqj and F idk = & 1 F i jdq j dq k ) 

The value of F that should be used for a particular physical problem therefore depends 
on what discretization convention one picks to use. Since the value of F is ambiguous, one 
approach might be to simply pick a discretization convention (it doesn't matter which one), 
and then choose F however necessary to reproduce the desired equilibrium distribution (ffl). 
To be concrete, let's choose Stratanovich convention. The Stratanovich Langevin equation 
(prp|) is well known to be equivalent to the following Fokker-Planck equation for the time 
evolution of the probability distribution P(q, t): 

P = V, [Te ia V 3 {e ja P) + F^P] . (14) 

My requirement is that this equation have the equilibrium distribution (^) as a time- 
independent solution. A simple way to achieve this is for 

Te ia Vj{e ja P e(l ) + F^P eq = 0, (15) 

which gives 



1 For a general review of background material for this paper, in notation close to what I use here, see, for 
example, chapter 4 of ref . Q| . The most substantial differences in notation are that my F and T are that 
reference's if and ifL 
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(11) 



F*** = (a-^VjV - Te ia e jaJ . (16) 
Compare to our naive starting point (§). Equivalently, 

F- t6 = {cr~ l )ijVjV — T(a~ 1 ) i jj . (17) 



So one might suspect that the Stratanovich equation (|10D together with (|16|), or equivalently 
the Ito equation @ together with (fT7|), gives the correct description of the system. As we 
shall see, this is indeed the case. 

Based on the presentation so far, the reader might be suspicious of two things. First, 
once we give ourselves license to go mucking around with the equation, changing what we 
thought was F(q) to suit our needs, how do we know we aren't supposed to change other 
things as well? In particular, what tells us that we shouldn't change e(q) in some way, then 
make some compensating change in F to force the equilibrium distribution to work out? The 
basic answer is that F is sensitive to the details of short-time physics and regularization, 
whereas e is not, as evidenced by the fact F must be changed when one adopts different 
conventions like Ito or Stratanovich, but e need not. This boils down to a discussion of 
the renormalizeability of the theory, and its consequences for how, in principle, the theory 
should be matched to a more fundamental description of the physics. Such issues are more 
familiar in the context of theories defined by path integrals, and so much of the rest of this 
paper will be to translate the discussion into that language. 



Secondly, the Stratanovich F of ([R]) is not the unique solution to ([14]) for P = P eq . The 
general solution is 

F st rat = (r-i^.y - Te m e jaj + h t e +y / T , (18) 

where h = h(q) is any function with V • h = 0. 

I shall throughout focus on systems where the underlying physics, whatever it may be, is 
reversible. More specifically, I shall assume that the effective theory ([!]) must be defined in 
such a way that equilibrium time-dependent correlations, such as (q(t) q(0)), are invariant 
under time reversal. [As discussed in Appendix [A|, the behavior of the Langevin equation 
with inertia (|5[) has this property, even though the equation by itself is not time-reversal 
invariant.] As I'll show, this assumption will rule out the extra term involving h in (|i~8"|). I 
shall then discuss in more detail, in the language of path integrals, how the issue of whether 
e or F should be modified from their "naive" values is an issue of renormalizeability. Specif- 
ically, I will define "naive" by assuming that there is some more fundamental description of 
the effective theory, such as the Langevin equation (|5|) with inertia, that is unambiguous and 
is described in terms of the same degrees of freedom q. The "naive" values of e and F will 
simply be those defined by naively taking the low frequency limit of the more fundamental 
equation. I'll show that terms associated with e require no ultraviolet renormalization, while 
terms associated with F do. Based on general procedures for matching effective theories to 
more fundamental underlying theories, one may then argue that F, and not e, should be 
modified to make the physics work out right. This is something that, I believe, may be well 
known to the few people to whom it is well known. However, since there seems to be general 
confusion on this matter, it seems worthwhile to continue rather than simply ending here. 
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III. THE PATH INTEGRAL VERSION 



The path integral corresponding to the Stratanovich Langevin equation (|T0|) is, somewhat 
imprecisely, 1 



7>(q",q',t"-^ 



q(i")=q" 
q(t')=q' 



[tZq(t)] exp 



(it L(q, q) 



(19) 



1 11 T 

£(q, q) = (q + F)i°ij(q + F)j - - F iyi + - e ka>k (q + F)i + — e ia>j e ja>i + 5(0) trine. 

(20) 

Here and henceforth, I abbreviate F Strat as simply F. The imprecision is just due to the 
fact that this path integral depends on the details of how time is discretized. I've implicitly 
assumed a time-symmetric discretization above. Specifically, (|T9|) really means 



/■q(t")=q" 
P(q",q',t" -t') = lim N \ 

^ -0 Ja(t')=Q' 



At- 



'q(*')=q' 
x exp 



I t 



-At$> 



qt - qt-At <lt + qt-At 
At ' 2 



(21) 



where iV is an overall normalization I shan't bother being explicit about, and where 5(0) in 
(ffiH) means (At) -1 . If I had used some other discretization convention in the path integral, 
the Lagrangian L would be different from (EDj) . 

For my argument, it will be sufficient (and convenient) to focus on fluctuations about 
equilibrium. It is then enough to consider the t' — > — oo and t" — > +oo limit of the path 
integral, writing 



[ctq(t)] exp 



dt L(q, q) 



(22) 



Because the system is dissipative, the boundary conditions on q at t = ±oo decouple. 
Equilibrium correlation functions like (q(t)q(O)) can then be evaluated in the usual way as 



(•■■) = Z- 1 [ [dq(t)\ exp [- / +0 ° dt L(q, q) 

J L J —CO 



(23) 



This form allows us to rewrite the Lagrangian in a form where the assumed time-reversal 
symmetry of correlations is manifest. Simply as an illustrative example, suppose for the 



2 There is a small change of notation from ref. 0: the (dete) 1 pre-exponential factor of that reference 
has been moved up into the exponent and absorbed into L. 

3 The formalism here is somewhat similar to that used to describe stochastic processes on curved manifolds, 
if one were to interpret cry as a metric tensor gij. I should emphasize that this is not the problem I 
am studying. On a curved manifold, the desired equilibrium distribution would be [det cr] 1 / 2 exp(— V/T) 
instead of exp(~V^/T). Mathematically, one can convert between these two problems by setting Vmanifoid = 
V + \T In deter. 
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moment that F did have the form fllTf ) that I'm trying to demonstrate. The Lagrangian 
( P0f) is not manifestly invariant under t — > — t. However, for the choice fli"6l) of F, this 
Lagrangian can be rewritten as 

L(q,q) = ^q l a lJ q J + d V) - V,V] + -(a 

+5(0) trine + (24) 

Only the last term is not manifestly time- reversal invariant. However, with the path integral 
in the current form of (0), we are allowed to throw away terms in L that are total time 
derivatives. That's because they can be integrated and re-expressed in terms of the boundary 
conditions, which we know are irrelevant. [As reviewed in Appendix [B], the fact that one 
may naively integrate total time derivatives is dependent upon the use of the symmetric 
time discretization ( pi)) .] The upshot is, that for the purpose of computing equilibrium (but 
time-dependent) correlators, we can replace L by 

£a(q, q) = ^ + ^ ( W)(a-%(V^) - iv^- 1 )^^] + ^(a l ) m 

+5(0) trine. (25) 

It's important to emphasize that L 2 cannot be used in the finite-time path integral (p!9[): it 
would not produce the same "P(q", q', t" — t'). 

In passing, we can now see the problem if F contained the additional h term of (p~8|) . 
The Lagrangian L given by (p0|) would then produce additional terms in fl24|) , one of which, 



— q i a ij h j e VIT , (26) 

is odd under time reversal. So the appearance of h would contradict my assumption about 
reversibility, unless fl2"6|) can be waived away as a total time derivative. This would require 

5Fi = h ie y/T = (a-%V j( j) (27) 

for some function 0(q). But, returning to (|T6|), this just corresponds to a shift of V — > V + (fi, 
and the equilibrium distribution produced by the Fokker-Planck equation ([TJ]) would then 
be exp[— {V + 4>)/T] instead of the required exp(— V/T). 

Now let's divorce ourselves from the particular form QIED that I've claimed for F. Instead, 
just write down the most general form for a manifestly time-reversal invariant L 2 : 

^(i^^ftMil^ + W ( 28 ) 

Before discussing how to determine the correct choices of and U for this equation, I 
want to discuss its renormalization.[] Imagine that the Lagrangian could be expanded in 



4 If the discussion in this paper is applied to a field theory, then I am assuming that the theory has already 
been regulated spatially, by the introduction of a small distance cut-off (either explicitly, by a spatial lattice 
cut-off, or implicitly by, for example, dimensional regularization) , and that I am now separately considering 
the issue of regularizing and renormalizing small-time divergences. The discussion of divergences and power 
counting will therefore be somewhat different from what is usual in field theory, where one typically regulates 
space and time simultaneously. 
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perturbation theory about some qo with cr (qo) finite and non-zero. For simplicity, say 
qo = 0. The arguments that follow will be about the entire perturbative expansion, to all 
orders in perturbation theory. To this end, write 

■£p ft <*ij (q) ft = 7^ [ft °ij (0) ft + ft Sa ij (q) ft'] > ( 29 ) 

and formally consider 5a SIS db perturbation. 

I will now show that the U(q) term in ( pffi ) requires counter-terms, but the c(q) term 
does not. The Lagrangian (^) is super-renormalizeable: counter-terms are needed only for 
one loop and two loop diagrams. This can be seen by simple power counting of diagrams. 
The most ultraviolet (UV) divergent diagrams will be those whose vertices involve as many 
powers as possible of internal frequency. That is, vertices generated by the ft Saij(q)qj terms 
in fl25|) will give more UV divergent behavior than vertices generated by the U(q) terms. 
And the most divergent behavior will occur when the q's are associated with internal rather 
than external lines. The most divergent one-loop diagrams are therefore those depicted by 
fig. |I|. The UV behavior of such graphs is 

J duj J^jxy = J du, (30) 

where {uj 2 ) v are the powers of u from the V vertices q<5crq, and {uj 2 )^ 1 are the uo — ► oo 
behavior of the I = V internal propagators. So this class of diagrams is linearly divergent in 
the UV. The linearly divergent piece is independent of external frequencies, and so requires 
a counterterm in U(q) but not one in the kinetic term qcrq. That is, so far we have seen 
that the effective U(q) is dependent on the details of short-time physics, but we have not 
seen any evidence that <r(q) is. 




FIG. 1. The UV divergent 1-loop graphs. 



Are there divergences associated with external frequency dependence? Expanding the 
diagrams of fig. [I] in powers of external momenta, each such power will come at the expense 
of one of the factors of u in (|3^) . So no counter-terms involving two external frequencies are 
necessary — that is, no counterterm for qcrq is necessary. There could be a log divergence 
associated with one external frequency factor, except that the resulting u integral / duo/uj 
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vanishes by the assumption that the equilibrium physics is time-reversal invariant (and the 
fact that I am using a symmetric discretization of the path integral, which respects that 
invariance) . 

Similarly, one could imagine a different method of hooking up the diagrams of fig. [I] in 
which one or more of the q factors were associated with external rather than internal lines. 
If this happens for two or more factors, the loop integral is no longer UV divergent and 
so in not sensitive to short-time physics. If it happens for only one factor, then the loop 
integral again vanishes by the assumption of time-reversal invariance. Therefore, the linear 
UV divergence ( |5D| ) previously discussed is the only UV divergence at one loop. 

Now consider two loop diagrams, again maximizing the divergence by using q 5o~ q ver- 
tices rather than U(q) ones, and associating all q's with internal lines. The result is shown 
in fig. |2|. Simple power counting gives the maximum possible new divergence at two loops 
as order 

/^w=/^' (3i) 

where I — V + 1 and to just counts powers of combinations of the two loop frequencies U\ 
and u)2- A logarithmic divergence is possible, but this divergence has no dependence on 
external frequencies. As before, if we focus on the dependence of these diagrams on external 
frequency, then there is no sensitivity to short-time physics. If we go to yet higher orders 
in diagrams, there are no new divergences at all: naive power counting gives J uj~ 2 ^ l ~^ d L cu 
for an L-loop diagram. 




\ 5a 



da' 

FIG. 2. The UV divergent 2-loop graphs. In the vertices with three internal lines, the q's can be 
associated with any two of them. 

The upshot of all this is that only U (q) is sensitive to the details of short-time physics; 
qaq is not. 

In the language of path integrals, the "naive" low frequency limit of a theory is what 
you get if you simply replace the Lagrangian by its low frequency limit. This replacement 
is naive when there are UV divergences in the resulting effective theory, which make the 
physics of that theory sensitive to how it is cut off in the UV — that is, sensitive to the 
details of the more fundamental theory. Fortunately, these details can be absorbed into a 
redefinition of the effective interactions which depend on them. In this case, that means 
U(q) is modified from its naive form, but qcxq is not. 
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Since this point is important, let me state it another way. The basic theory and technol- 
ogy of matching parameters of low-frequency effective theories to those of more fundamental 
theories, to any and all orders in perturbation theory, has been used in a number of problems 
in field theory, including Bose condensation Q, ultrarelativistic plasmas ||, heavy quark 
physics 0, and non-relativistic plasma physics 0. (For a general discussion, see also ref. 
||.) The basic idea is to fix a regularization for the effective theory, to then treat all the 
parameters of the effective theory as adjustable, to calculate low-frequency observables in 
both the effective theory and the more fundamental theory, and then to fix the parameters 
of the effective theory to obtain agreement. If we simply fixed cr(q) and U(q) to be their 
naive values, everything would match order by order in perturbation theory if it weren't 
for the UV divergences in the effective theory — if diagrams only involved low frequencies in 
propagators and vertices, then each internal line and vertex in the underlying theory would 
match up with each one in the naive effective theory. The only thing we have to correct 
for are UV divergent graphs or sub-graphs, where the correspondence would fail. But the 
previous discussion shows that such divergences have the form of counter-terms for U(q), 
and so the matching of the two theories can be fixed by appropriate adjustment of those 
counter-terms. 

If we had a specific underlying theory in mind, we cold carry out this matching procedure 
by computing to two loops in perturbation theory. The crucial point here is that such a 
calculation is unnecessary, and detailed knowledge of the underlying theory is unnecessary, 
if we know that the equilibrium distribution is (f|). We can then immediately solve the 
problem by simply requiring that we choose U(q) so that the equilibrium distribution comes 
out right. 

Returning to the Lagrangian fl20D corresponding to the original Langevin equation, the 
only way we can change U(q) without changing qcrq is by changing F, which justifies the 
argument for choosing F Strat given in section J|. Alternatively, one can uniquely determine 
U(q) directly by requiring that the Euclidean Schrodinger equation, corresponding to the 
path integral with the generic Lagrangian L 2 (p8[) , generate exp(— V/T) as its equilibrium 
distribution. This is carried out in Appendix |C|, which also explains some interesting dis- 
tinctions between the interpretation of the Euclidean Schrodinger equation associated with 
L 2 [ fl25|) or (f28|) l and that associated with L (20). The final result agrees with the initial 
analysis of section (|II|). 
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APPENDIX A: THE LANGEVIN EQUATION WITH INERTIA 

The Langevin equation (|5|) with inertia can be written in the form of a standard first- 
order Langevin equation by rewriting it as a system of first order equations, introducing 
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v + <x(q)v = -W(q) + e- 1 (q)£, (Ala) 
q = v. (Alb) 

This equation is free of the usual discretization ambiguities. To see this, consider the differ- 
ence between Stratanovich and Ito discretization conventions, which differ in how the q in 
the e -1 (q)£ term is evaluated. That difference in q is |qAt, which means a difference in 
the e _1 (q) £ term of 



At ( ^_ x 
~2 



e v )ij,k(lkij- (A2) 



In a Langevin equation such as our effective theory (|6]), this would not vanish in the At — ► 
limit because £ is order (At) _1//2 , and then so is q by the equation of motion (|6p.Q In the 
case at hand, matters are different, because q k = v k , and, though ( |Ala|) implies that the 



time derivative v will be order (At) -1 / 2 , v itself is perfectly finite as At — ► 0. So, in the 
case at hand, the difference ( |A~^ ) vanishes in the At — > limit. 

It is also interesting to sketch how this works in the path integral formulation. Converting 
equation ([5]) directly into a path integral using standard methods[] gives a path integral with 
Lagrangian 



L = — 



(q + aq + W)<r _1 (q + aq + W) + L v + 5(0) tr In e, (A3) 



L n = rji {dijdf + Oijdt + (e x ) ai e kaJ q k + (e l ) a k,Ak + ^i{e ka S7 k V) (A4) 

where fj and 77 are anti-commuting Grassman fields (ghosts). A diagrammatic analysis 
similar to that discussed in the main text produces no UV divergences because the qcr _1 q 
term means that q propagators go like u>~ 4 in the UV instead of only lu~ 2 . So diagrams 
converge faster in the UV than they did in the previous analysis. 

Now let's check the equilibrium distribution for q. Since (|A1|) formally has the same 
form as a generic first-order Langevin equation (||) [the number of degrees of freedom have 
simply doubled from q to (q, v)], we can use the result flTil) for the Fokker-Planck equation 
for the probability distribution P(q, v,t), which translates in the present case to 



P = V Vi {T(e- 1 ) ja V t)j [(e- 1 ), a P] + ((TV + V q V),P} + V ft {-v.P} 

rv v y V] p + v Vi ■ (vjP)] + v q y • v v p - v ■ v q p. 



The equilibrium solution is 



(A5) 



5 Interestingly, this is a very special case of a general set of Langevin equations considered in ref. 

6 See, for example, sees. 4.7-8 of ref. ^ for a review. 

7 See, for example, sec. 4.6 of ref. Q for a review, as well as ref. ||. 
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P cq (q, q ) = exp {-1 [\v 2 + V(q)} } . (A6) 

If we are only interested in the distribution in q and not v = q, we may integrate over all 
possible values of q and obtain (f|). 

Finally, I turn to my claim that this system has time-reversal invariant equilibrium 
correlations (q(t)q(O)). One way to see this is from the path integral Lagrangian L of (|A3|) , 
which can be rewritten as 

L = h { ^ + vv)a ~ l {il +vv) + h ^ +Lr > + m tr ln e + h dt + v ) ■ (A7) 

As discussed in the main text, the total time derivative at the end can be thrown away if 
one is only interested in equilibrium correlations. The rest of L is manifestly time-reversal 
invariant, except possibly the ghost term L v . (It is again important here that my path 
integrals are all implicitly defined with symmetric time discretization.) 

For this system, the ghost term is actually fairly trivial, as can be seen by integrating 
out the ghosts. Consider the perturbative expansion of the theory, and imagine an arbitrary 
ghost loop. In frequency u, the pertrubative ghost propagator is of the form [— u 2 — ia(0)uj + 
M 2 ]~ l (where I have suppressed the indices The poles in uo all lie in the lower half 

complex plane. Almost any frequency integral corresponding to a ghost loop can then be 
seen to vanish by closing the integration contour in the upper half plane. The one exception 
is for loops with frequency integrals where the integrand does not fall off faster than 1/ou as 
\u\ — > oo, since in this case we cannot ignore the contribution of the contour at infinity. This 
situation arises only for ghost loops with a single ghost propagator and a single vertex of 
the form — iu5a(q). The resulting frequency integral yields an effective interaction among 
the q's of 

r+°° du -iuSa x 

Loo 2^ [-u* - iu a(0) + M*] ~ ~ 2 { ' 

Restoring indices i,j and not worrying about keeping track of additive constants in the 
action, this becomes, 

L v -> -§tr<7. (A9) 
This is clearly time-reversal invariant. 



APPENDIX B: INTEGRATING BY PARTS WITH SYMMETRIC 

DISCRETIZATION 



Consider a term in the action of the form 

/ = J dtq.V/(q), (Bl) 

which can be naively integrated to yield a boundary term. I will review how this naive 
integration is justified for the symmetric time discretization used in this paper. In that 
discretization, I really represents 
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/ = £(q. + i-q.)-V/(^-^y (B2) 



where r = tj At is an integer parameterizing the time steps. Now consider the Taylor 
expansion of /(q r +i) — /(<1t) about q = (q r +i + qr)/2: 

/(qr+i) - fM = Aq • V/(q) + 0[(Aq)% (B3) 
where Aq = q r +i — q r . So we can rewrite 

* = E {/(Qr+i) - /(qr) + 0[(Aqf]} . (B4) 

r 

In a path integral, contributions to the action survive in the At — > limit if they contribute 
O(At) or more per time step.| The kinetic term determines the size of Aq to be Of(At) 1 / 2 ]. 
Therefore, the 0[(Ag) 3 ] term in (|B4[) can safely be ignored, and what's left trivially cancels 



between successive time steps, except for boundary terms. 

The fact that the error was 0[(Ag) 3 ] and not 0[(Aq) 2 ] (which would not be ignorable) 
in (]B4]) depended crucially on the symmetric discretization q = (q T+ i + q T )/2. 

APPENDIX C: EUCLIDEAN SCHRODINGER EQUATION FOR L 2 

We can directly determine the equilibrium distribution generated by the generic La- 
grangian (|28|) by transforming the path integral into a Euclidean Schrodinger equation. 
First, recast the path integral over q(t) as a path integral over q(t) and momentum p(t): 



Z 2 = lim N' 

At-,0 



11 dp T dq T 



5a(p,q) = E {-<Pr ■ (q T " Qr-i) + AtF 2 (p T , Sl±^i^| , ( C2 ) 



^ 2 (p,q)=Tp ( T- 1 (q)p + M (q), (C3) 

u = U- ^trlne, (C4) 

which can be checked simply by doing the Gaussian integration over p. In this form, the 
path integral is well known to correspond to a Schrodinger equation 

[#2(p,q)W>(q,t) = -^(q,*), (C5) 



8 As a quick mnemonic, think about the contributions of a potential term: J dt U(q) = At J^. C/(q). Such 
contributions are manifestly O(At) per time step. 
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where the subscript W indicates Weyl ordering of the operators p and q. The Weyl ordering 
formula we need in the case at hand is that[] 

\p iPj A(q)] w = ^ {pi, {p v A{q)}} , (C6) 

where I've now introduced hats to emphasize that p and q are operators. Taking p = — iV, 
the Schrodinger equation ( Pop then becomes 



i> = \ t 



[<y )ijVjVj + (a )ij,{Vj + -{a )ij/n 



1 

4 l 



uU. (C7) 



I have used the symbol ^ instead of P in this equation because the equilibrium result for 
ip must be the square root of the probability distribution. That is, we want the Schrodinger 
equation ( |C7| ) to have 

^ eq = e~ v ^ T (C8) 

as its solution rather than exp(— V/T). To understand this, consider the original Lagrangian 
L from the discussion of the path integral form of Langevin equations. L differed from L 2 
only by a total time derivative. The Schrodinger equation corresponding to L is simply the 
Fokker-Planck equation flHj), which I'll now write in the form 

P = -H X P. (C9) 

Both formulations, in terms of L or L 2 , should generate the same equilibrium physics — that 
is, the same time-dependent correlation functions. A crucial difference between H 2 and Hi, 
however, is that the operator H 2 is Hermitian, while Hi is not. Now think about what it 
means to write down a path integral expression for the equilibrium probability distribution 
P eq (q) = (5(q — q)) in terms of actions that run between arbitrarily large times — T and 
+T. The corresponding object in the Schrodinger formulation is 

p (a) = hm (q(+^)l^ r lq)(ql^ T |g(-^)) (C10) 

eqV ; r-°o ( q ( +T )|e- 2 ^|q(-T)) ' V ; 

where H can be either Hi or H 2 . What dominates the long-time evolution operator 
exp(— HT) is the equilibrium state, which I'll denote |eq). The difference between Hi and 
H 2 is that the long-time evolution generated by the H 2 must be symmetric in its overlap 
with the initial and final states, becuase H 2 is Hermitian. That is, in the large T limit, 

(q'|e-^ r |q // )^(q / |eq 2 )(eq 2 |q"). (Cll) 

Hi is not Hermitian and so does not have this symmetry. In fact, we know from the usual 
Fokker-Planck equation corresponding to Hi that the evolution is dissipative, and the result 
of long-time evolution is independent of initial conditions. So 



See, for example, ref. || for a review of this fact in the present context. 
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(q'|e-^|q"> -> (q'K). 
In the case of H 2 , (|U10 ) then becomes 

2 

^eq(q) = (q|eq 2 ) , 

whereas for Hi it becomes 

A q (q) = (q|eq x ). 



(C12) 



(C13) 



(C14) 



One can now see that the equilibrium wave function represents the square root of P eq in the 
case of H 2 but P eq itself in the case of H 1 . 

In any case, we can now uniquely determine w(q), and hence U(q), simply by requiring 
that the equilibrium amplitude ( |C8D be a solution of the Schrodinger equation (p7|). One 
finds 



u 



^(VV)a-\VV) - \v{a- l VV) + 



(C15) 



which precisely reproduces the result ( p5|) for the L 2 that describes the Langevin equation, 
with my claimed result (fl6|) for F strat . 
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